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The application of the quasi-steady-state approximation to the Michaelis-Menten reaction embed- 
ded in large open chemical reaction networks is a popular model reduction technique in deterministic 
and stochastic simulations of biochemical reactions inside cells. It is frequently assumed that the 
predictions of the reduced master equations obtained using the stochastic quasi-steady-state ap- 
proach are in very good agreement with the predictions of the full master equations, provided the 
conditions for the validity of the deterministic quasi-steady-state approximation are fulfilled. We 
here use the linear-noise approximation to show that this assumption is not generally justified for 
the Michaelis-Menten reaction with substrate input, the simplest example of an open embedded 
enzyme reaction. The reduced master equation approach is found to considerably overestimate the 
size of intrinsic noise at low copy numbers of molecules. A simple formula is obtained for the rela- 
tive error between the predictions of the reduced and full master equations for the variance of the 
substrate concentration fluctuations. The maximum error is reached when modeling moderately or 
highly efficient enzymes, in which case the error is approximately 30%. The theoretical predictions 
are validated by stochastic simulations using experimental parameter values for enzymes involved 
in proteolysis, gluconeogenesis and fermentation. 



It is well known that whenever transients in the con- 
centration of a substrate species decay over a much 
slower timescale than those of the enzyme species, one 
can invoke the quasi-steady-state approximation (QSSA) 
to considerably simplify the deterministic (macroscopic) 
rate equations. 1,2 The study by Rao and Arkin 3 pio- 
neered the use of the same approximation but on a meso- 
scopic level, i.e., applying a stochastic version of the 
approximation to obtain reduced chemical master equa- 
tions. This approximation has since become ubiquitous 
in stochastic simulations of large biochemical reaction 
networks inside cells (see, for example, Refs. 4-7) al- 
though its range of validity is presently unknown. A 
plausible hypothesis is that the stochastic QSSA is valid 
in the same regions of parameter space where the deter- 
ministic QSSA is known to be valid. A handful of numer- 
ical studies 8,9 have shown that for some choices of rate 
constants which are compatible with the deterministic 
QSSA, the differences between the reduced and full mas- 
ter equation approaches are practically negligible. How- 
ever, none of these studies exclude the possibility that 
there exist regions of parameter space where the deter- 
ministic QSSA is valid but the stochastic QSSA exhibits 
large systematic errors in its predictions. In particular, 
one is interested in knowing how accurate are the predic- 
tions of the stochastic QSSA for the size of intrinsic noise, 
i.e., the size of fluctuations in concentrations, since such 
noise is known to play important functional roles in bio- 
chemical circuits. 10 Numerical approaches cannot easily 
answer such questions because the stochastic simulation 
algorithm, the standard method which exactly samples 
the trajectories of master equations, 11 is computationally 
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expensive. 

In this communication, we seek to develop a theoret- 
ical approach to answer the following question: Given 
that the rate constants are chosen such that the deter- 
ministic QSSA is valid, what are the differences between 
the predictions of the reduced and full master equations 
for the variance of the fluctuations about the mean con- 
centrations? We obtain a formula estimating the size 
of these differences for the simplest biochemical circuit 
which embeds the Michaelis-Menten reaction and con- 
firm its accuracy using stochastic simulations. We find, 
using physiological parameter values, that the reduced 
master equation approach can overestimate the variance 
of the fluctuations by as much as ~ 30%. 

We start by considering the Michaelis-Menten reaction 
with substrate input 

Xs, Xs + Xe Xc Xe + Xp, (1) 

ki 

where Xj denotes chemical species i and the fc's denote 
the associated macroscopic rate constants. The reaction 
can be described as follows. Substrate molecules (species 
S) are pumped into some compartment at a constant 
rate, they bind to free enzyme molecules (species E) to 
form substrate-enzyme complexes (species C) which then 
either decay back to the original substrate and free en- 
zyme molecules or else decay into free enzyme and prod- 
uct molecules (species P). The first reaction in (1) could 
equally represent the production of substrate by a first- 
order chemical reaction provided the species transform- 
ing into substrate exists in concentrations large enough 
such that fluctuations in its concentration can be ignored. 
The sum of the concentrations of free enzyme and com- 
plex is a constant since the enzyme can only be in one of 
these two forms. Hence, all mathematical descriptions of 
the Michaelis-Menten reaction can be expressed in terms 
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of just complex and substrate variables. On the macro- 
scopic level, the QSSA proceeds by considering the case 
in which transients in the complex concentration decay 
much faster than those of the substrate. This condition 
of timescale separation is imposed by setting the time 
derivative of the macroscopic complex concentration to 
zero, solving for the steady-state complex concentration 
and substituting the latter into the rate equation for the 
substrate concentration which leads to the new rate equa- 
tion 

d r y m] , k 2 [E T )[X s {t)} 

=**"- K M + [Xs(t)] ' (2) 

where [Xg(t)] is the substrate concentration at time t, 
Km = (k\ + 1^2) /ko is the Michaelis-Menten constant, 
and [Et] is the total enzyme concentration, i.e., the sum 
of the concentration of free enzyme, [Xe(£)], and of the 
concentration of complex, [Xc(t)], which is a constant 
as previously mentioned. Note that the notation [Xj\ 
(without explicit dependence on t) will be reserved for 
the steady-state concentration of species i. We assume 
that at t = wc are in steady-state conditions; note that 
the system (1) is guaranteed to have a stable steady- 
state if the condition km/k^Ex] < 1 is satisfied. Linear 
stability analysis of the full rate equations describing (1) 
shows that the timescale for the decay of transients in the 
substrate concentrations is t s = (Jco[Xe]) , while the 
timescale for the decay of transients in the complex con- 
centrations is t c = (k [Xg]+ki + k 2 )^ 1 ■ Hence, the crite- 
rion for the validity of the QSSA on the macroscopic rate 
equations (the deterministic QSSA), i.e., for the validity 
of Eq. (2), reads r s /r c = 7 = ([X S }+K M )/[X E ] > 1 (see 
also Ref. 13). This condition implies that timescale sep- 
aration is possible over the whole range of macroscopic 
substrate concentrations whenever the total enzyme con- 
centration is much smaller than the Michaelis-Menten 
constant. 

The question that we address in the rest of this com- 
munication is the following: Given that the condition 
7 ^> 1 is satisfied, what is the variance of the noise about 
the mean concentrations as predicted by the reduced and 
full master equations? 

The stochastic QSSA method implicitly starts by de- 
ducing that Eq. (2) is effectively the rate equation one 
would associate with a system of two chemical processes 

^ Xg, X s ^ X P , (3) 

where k' is an effective (time-dependent) rate constant 
equal to k 2 [E T ]/(K M + [Xs(t)]). Note that while the first 
reaction is elementary, the second is clearly not, since it 
can clearly be broken down into a set of more funda- 
mental constituent reactions. Given the reduced set of 
reactions (3) one can then construct a reduced master 
equation for the set of reactions (1) (see Ref. 14 and 
supplementary material for the construction of master 



equations) 

n 

-P(n s ,t) =nki n (E% 1 - l)P(n s ,t) 

+ {Er-l)^^--P{n s ,tl (4) 
KM + ns/vl 

where f2 is the compartment volume in which the reac- 
tions are occurring, ng is the absolute number of sub- 
strate molecules, P(ng, t) is the probability that the sys- 
tem has 715 substrate molecules at time t and E™ is 
the step operator which upon acting on a function of 
ng changes it into a function of ng + m (Ref. 14). Wc 
note and emphasize that the physical basis of this master 
equation is not clear because such equations have been 
derived from first principles for elementary reactions, 16 ' 1 
while (3) involves a non-elementary reaction. Equa- 
tion (4) is simply written by analogy to what one would 
write down for (3) if both reactions were elementary and 
hence its legitimacy is a priori doubtful. 

Now we want to use this master equation to deduce 
the variance of the noise in the macroscopic substrate 
concentrations. It is well known that in the macroscopic 
limit, the master equation for monostable chemical sys- 
tems can be approximated by a linear Langevin equa- 
tion, an approximation called the linear noise approxi- 
mation (LNA). 14 ' 19 For systems with absorbing states or 
exhibiting multimodality, the LNA will not usually give 
accurate results (see, for example, Ref. 18) but its ap- 
plication to our example, the Michaelis-Menten reaction 
with substrate input, is not problematic since this reac- 
tion is only capable of monostable behavior. The steps 
to construct the LNA for a general monostable chemical 
reaction system are summarized in the supplementary 
material. Here, we will simply state the results of this 
recipe when applied to the master equation, Eq. (4). The 
Langevin equation approximating Eq. (4) in the macro- 
scopic limit of large molecule numbers and in steady-state 
conditions is 

m ns{t) = --ns® + + -j-)ntl (5) 

where ng (t) denotes the fluctuations about the macro- 
scopic steady-state substrate concentration defined as 
T)s(t) = ng{t)/Q, — [Xs] and T(t) is white Gaussian noise 
defined by (T(t)) = and (T(t)T(t')) = S(t - f). From 
Eq. (5), one can show (see supplementary material) that 
the variance of concentration fluctuations 775 (t) in steady- 
state conditions is given by 

[Xg]f [X S ]\ 

where the subscript "sLNA" stands for "LNA of the mas- 
ter equation reduced using the stochastic QSSA". Note 
that [Xg] in Eqs. (5) and (6) is the steady-state substrate 
concentration obtained by solving for [Xg] from Eq. (2) 
with time derivative set equal to zero. 
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We shall now derive expressions for the variance of 
substrate concentration fluctuations using the full mas- 
ter equation approach. The steps of this method are as 
follows: (i) one writes down the master equation for the 
elementary chemical processes (1), (ii) the two Langevin 
equations for the complex and substrate fluctuations are 
obtained using the LNA of the master equation, (iii) ex- 
pressions are found for the variance of complex and sub- 
strate concentration fluctuations in steady-state condi- 
tions, (iv) the limit 7 > 1 is taken of the expressions 
derived in step (iii), leading to the final expressions for 
the variance of substrate fluctuations about the steady- 
state substrate concentration solution of the rate equa- 
tion, Eq. (2). We note that this method, unlike the first 
one, does not make any assumptions about the validity 
of an ad-hoc reduced master equation since it is based on 
the master equation for elementary processes and hence 
is guaranteed to be correct. 20 We now proceed to put this 
systematic recipe in practice. 

The master equation for the four elementary chemical 
processes given by (1) is 



d t P(n s ,n c ,t) 



^(E^Ec 1 -l)n s (n T -n c ) 
VLk m {Eg X - 1) + k 1 (Eg 1 E+ 1 - l)n c 



l)n c 



P(n s ,nc,t), 



(7) 



where nc is the absolute number of complex molecules 
and tit is the absolute total number of molecules of en- 
zyme in free and complex form. Note that nr is a con- 
stant equal to [Er]^- In the macroscopic limit, the mas- 
ter equation, Eq. (7), can be approximated by a pair of 
Langevin equations as given by the LNA (see supplemen- 
tary information) 



7ft 



vc(t) 
vs(t) 



-{K M 
K X 4 



f- [Xs]) 
[Xs] 



[Xe] 



vc{t) 
vs{t) 



V k Q [X E ][X s ] - y/fci[X c ] 
fk~ -y/k [X E ][X s ] ,/k^Xc] 



-VHXc\ 




r(t), 

(8) 



where K\ = fei/fco, Vc{t), and rjs{t), respectively, denote 
the fluctuations about the macroscopic steady-state com- 
plex and substrate concentrations and T(t) is a 4 dimen- 
sional vector whose entries are white Gaussian noise with 
the properties = and (r^i)^-^')) = Si,jd(t-t'). 

Using Eq. (8), one can show (see supplementary informa- 
tion) that the variance of substrate concentration fluctu- 
ations rjs{t) in steady-state conditions is given by 



&LNA 



[Xs 



n 



7>1 



[X S ] K x + [X s ] 7 
K M Km + [X s ] 1 + 7 
1 , [X s ] K x + [X s ] 



Km K m 



[Xs 



(9) 



where in the last step we took the limit of 7 ^> 1, cor- 
responding to the condition in which the deterministic 
QSSA Eq. (2) is valid. 



Comparing Eqs. (6) and (9), we see that the two are 
not generally equal to each other except in the case 
/3 = fca/fci -C 1. From a fundamental point of view, 
this disagreement implies that the reduced master equa- 
tion does not obey the generalized fluctuation-dissipation 
theorem of nonequilibrium physics 14 ' 21 and that hence it 
is flawed. More importantly, we observe that the con- 
dition p = k^/ki < 1 is not equivalent to the quasi- 
stcady-state condition 7 > 1. The former condition is 
consistent with the enzyme-substrate complex being in 
thermodynamic equilibrium with free enzyme and sub- 
strate, a condition which is difficult to uphold in open 
systems since they are characterized by nonequilibrium 
steady states. While the quasi-steady-state condition can 
easily be satisfied in open systems since it is only required 
that the total enzyme concentration is much less than the 
Michaclis-Mcnten constant. Hence, we can conclude that 
for open systems, the stochastic QSSA based on Eqs. (3) 
and (4) is not the legitimate stochastic equivalent of the 
deterministic QSSA. 

There are two possible hypothetical scenarios which 
would imply that the stochastic QSSA is perhaps still 
a very good general method to estimate the size of the 
concentration fluctuations. The first case would be if 
experimental evidence showed that for many enzymes it 
just happens that (3 < 1. The second case would be 
if experimental evidence showed no such restriction on 
/3 but nevertheless the difference between the variance 
prediction of the reduced and full master equations is so 
small as to be negligible. We now consider each case. 

A perusal of the experimental data available in the 
literature shows that there are very few studies which 
simultaneously report values of k\ and k%, the data re- 
quired to estimate (3. The vast majority of studies report 
values for Km, a considerable number report &2 and a 
small percentage report both fc 2 and Km- 22 Now the ra- 
tio k^jKM, frequently called the enzyme efficiency, 23 is 
defined as 



e 



fc 2 

K M 



k . 



(10) 



The recent study by Bar-Even et al. 22 based on mining 
the Brenda 25 and KEGG databases 26 concluded that for 
most enzymes 8 lies in the range 10 3 — 10 6 M _1 s _1 . It is 
also known that the association constant ko takes values 
in the range 23 10 6 — 10 9 M _1 s _1 . We can conclude from 
these two pieces of data and using Eq. (10) that the range 
of P for most enzymes is between and some number 
which is much greater than 1 and that hence on the basis 
of experimental data one cannot argue for the general 
validity of the stochastic QSSA. 

Of course, as previously mentioned, it could still hap- 
pen that even though there is no restriction on /3, that 
the variance as predicted by the stochastic QSSA and the 
true variance are negligibly small. We can test this hy- 
pothesis quantitatively by using Eqs. (6) and (9) to derive 
the fractional relative error e in the variance prediction 
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FIG. 1: (a) Plot of the fractional relative error, e, in the variance of fluctuations as predicted by the reduced master equation 
obtained from the stochastic QSSA versus the enzyme saturation parameter a. The solid lines show the theoretical prediction 
Eq. (11) for three different values of j3: 1 (blue), 2.8 (magenta) and 10 (yellow-green). These three cases are respectively 
consistent with the enzymes being Chymotrypsin or Malate dehydrogenase, Lactate dehydrogenase and a highly efficient enzyme 
such as /3-Lactamase. The data points show stochastic simulations using the Gillespie algorithm for reactions involving a total 
of 100 enzyme molecules and total enzyme concentrations in the nano and millimolar range (see supplementary information 
for details), (b) Plot of the fractional relative error e versus the total number of enzyme molecules tit for the case a — 0.5 
and /3 = 10. The data points are obtained from stochastic simulations. The solid line is simply a guide to the eye. The total 
number of enzymes is varied at constant total enzyme concentration. The ratio of substrate and complex decay timescales, 
7, is greater than 10 in all cases shown in (a) and (b) which enforces the validity of the deterministic QSSA. Details are in 
supplementary information. 



of the stochastic QSSA 

e _ Qlna - QsLNA _ -(1 - a)a(3 

<7 LNA ! + /?(! -a(l -a))' 

where a = fejn/(&2[-^r])j a non-dimensional quantity 
which can take values between and 1 as previously men- 
tioned in the discussion after Eq. (2). Furthermore, it 
can be shown using Eq. (2) that at steady-state one has 
[X E ] = [E T ](l — a), and [X c ] = [E T ]a, from which we 
can deduce that a is a measure of how saturated is the 
enzyme with substrate. Note that Eq. (11) shows that 
the relative error tends to zero as a — > and a — > 1 
and that hence the reduced master equation provides 
a correct prediction of the size of the substrate fluc- 
tuations whenever the free enzyme or complex concen- 
trations are very small (similar results have been ob- 
tained by Mastny et al. 2A for the Michaelis-Menten re- 
action with no substrate input; however their results are 
not for general a and /3 and do not enforce the valid- 
ity of the deterministic QSSA; see later for discussion). 
In Fig. 1(a), the solid lines illustrate the predictions of 
Eq. (11) for three different values of (3: (i) 1, (ii) 2.8, and 
(iii) 10. Case (i) utilizes experimental data for the en- 
zymes Chymotrypsin and Malate dehydrogenase with re- 
spective substrates Acetyl-L-tryptophan and NADH, 27,28 
while case (ii) is based on data for the enzyme Lactate 
dehydrogenase with substrate pyruvate. 29 These enzymes 
are respectively involved in proteolysis, gluconeogenesis 
and the conversion of pyruvate (the final product of gly- 
colysis) to lactate in anaerobic conditions. Case (iii) 
showcases the largest possible error made by the stochas- 
tic QSSA; this is consistent with a highly efficient en- 
zyme such as /3-Lactamase for which is of the same 
order of magnitude as the maximum possible association 
rate constant 23 k ~ 10 8 - 10 9 .s _1 A/ _1 . The theoretical 



predictions of our LNA based method are confirmed by 
stochastic simulations of the master equations, Eq. (4) 
and Eq. (7), using Gillespie's algorithm 11 [data points in 
Fig. 1(a)]. Note that the maximum possible percentage 
error is about 30%, which is significant. Also note that 
the maximum error in all cases is reached at a = 1/2, 
namely, when the enzyme is half saturated with sub- 
strate which occurs when the substrate concentrations 
are equal to the Michaelis-Menten constant Km (this is 
the case for most enzymes of the glycolytic pathway 31 ); 
for substrate concentrations much smaller or larger than 
Km, the error is negligible. 

The LNA is, strictly speaking, valid for large vol- 
umes or, equivalently, in the limit of large number of 
molecules 14 ' 30 and hence one could argue that our the- 
oretical formula Eq. (11) is of limited validity inside 
cells, where molecule numbers can be quite small. 32 Fig- 
ure 1(b) shows the results of stochastic simulations for 
the case a = 0.5 and f3 = 10 using a total number of en- 
zyme molecules tit varying between 1 and 100 molecules. 
Note that the error e is practically constant at 30%, the 
value predicted by the LNA and shown in Fig. 1(a). This 
suggests that the estimates provided by our method are 
accurate even for low copy number conditions. 

Our study has focused on the most common type 
of stochastic QSSA in the literature which is heuristic 
in nature and hence the question regarding its valid- 
ity. There are a class of alternative model reduction 
techniques 24 ' 33 based on singular-perturbation analysis 
(sQSPA and sQSPA-Sl) which are rigorous and whose 
validity is not under question. For the Michaelis-Menten 
reaction without substrate input, these methods lead to a 
reduced master equation of the same form as the heuristic 
stochastic QSSA whenever the free enzyme or complex 
concentrations are very small (see Table II of Ref. 24). 
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This implies that for such conditions the error in the pre- 
dictions of the stochastic QSSA should be zero, a result 
which is also reproduced by our method. However, note 
that though these concentration conditions can be com- 
patible with the deterministic QSSA they are not syn- 
onymous with it. The sQSPA methods do not lead to a 
reduced master equation for parameters consistent with 
the deterministic QSSA and hence cannot make state- 
ments regarding the accuracy of the heuristic stochastic 
QSSA in such conditions. Our contribution fills this im- 
portant gap by deriving an explicit formula for the error 
in the predictions of the stochastic QSSA, i.e., Eq. (11), 
for all parameters values consistent with the determin- 



istic QSSA. We finish by noting that a recent study by 
Gonze et al. s also studied the reaction system (1) using 
numerical simulations and found little difference between 
the predictions of the stochastic QSSA and the full mas- 
ter equation. The study used values of ft = 0.1 (see Table 
7.2 in Ref. 8) and hence in the light of our results, it is 
clear why they observed high accuracy of the stochastic 
QSSA. However, as we have shown, this is not the general 
case: many enzymes have large ft and hence discrepancies 
of the order of few tens of percent between the predictions 
of the reduced and full approaches will be visible when- 
ever substrate concentrations are approximately equal to 
the Michaelis-Menten constant. 
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Supplementary information 

A. General formulation of master equations for 
elementary and non-elementary processes 

Consider a general chemical system confined in a com- 
partment of volume Q and consisting of a number N of 
distinct chemical species interacting via R chemical re- 
actions of the type 

k ■ 

sijXi + . . . + s Nj X N T\jX\ + . . . + r Nj X N . (12) 

Here, j is an index running from 1 to R, X^ denotes 
chemical species i, s^j and are the stoichiometric co- 
efficients, and kj is the macroscopic rate of reaction. Note 
that these reactions are not necessarily elementary (uni- 
molecular or bimolecular reactions). If the jth reaction 
is elementary then its rate kj is a constant while if it is 
non-elementary kj is a function of macroscopic concen- 
trations. The general form of the master equation for 
both cases is 1 

Ii N 

<>WlI /; - ' - lU(n,0)P(rM), (13) 

where P(n, t) is the probability that the system is in a 
particular mesoscopic state n = (ni, un) t and n 2 ; is 
the number of molecules of the ith species. Note that 
Ef is a step operator - when it acts on some function of 
the absolute number of molecules, it gives back the same 
function but with rij replaced by rii + x. The chemical re- 
action details are encapsulated in the stoichiometric ma- 
trix Sij = r,j — Sij and in the microscopic rate functions 
fj(n, £^)- The probability that the jth reaction occurs in 
the time interval [t, t + dt) is given by Clfj(H, Cl)dt. 

For elementary reactions, the microscopic rate function 
takes one of four different forms, depending on the order 
of the jth reaction: (i) a zeroth-order reaction by which a 
species is input into a compartment gives fj (ft, fi) = kj ; 
(ii) a first-order unimolecular reaction involving the de- 
cay of some species h gives fj(ri,fl) = fcj-n/jfi -1 ; (hi) a 
second-order bimolecular reaction between two molecules 
of the same species h gives fj(n, 0) = kjrih(nh — l)^ -2 ; 
(hi) a second-order bimolecular reaction between two 
molecules of different species, h and v, gives fj(n, Q) = 
kjUhfiyil" 2 . Note that these forms for the microscopic 
rate functions have been rigorously derived from micro- 
scopic physics 2,3 and hence the validity of Eq. (13) for 
elementary reactions is guaranteed. 

For non-elementary reactions, the form of the mi- 
croscopic rate function has to be basically guessed by 
analogy with the prescription for elementary reactions. 
For example, for the set of reactions (3) in the main 
text, the second reaction is a non-elementary first-order 
reaction with a time-dependent macroscopic rate con- 
stant k'{t) = k 2 [E T ]/(K M + [X s (t)]), where [E T ] is the 
constant macroscopic total enzyme concentration and 



[Xs(t)} is the instantaneous macroscopic concentration 
of species S. Hence, one would use the microscopic rate 
function /(n, fl) = k 2 [ET](ns /£>)/ (Km + ng/Cl) based 
on the formula stated above for an elementary first-order 
reaction. Of course, master equations based on micro- 
scopic rate functions obtained from this procedure are 
ad-hoc and have no fundamental basis. 



B. General formulation of the linear noise 
approximation in steady-state conditions 

Here, we provide a step by step recipe to construct the 
linear noise approximation (LNA) of the master equa- 
tion, Eq. (13), for the set of reactions (12). We note 
that this approximation is only valid for a monostable 
system [the condition is formally given by Eq. (3.4) in 
Ch. X of the book by van Kampcn 1 ]. Let the macro- 
scopic steady-state concentration of species i be given by 
[Xj] and the derivative with respect to this variable be 
denoted by V^. Furthermore, we shall distinguish matri- 
ces by underlining them. The five steps to constructing 
the LNA in steady-state conditions (for both elementary 
and non-elementary reactions) are then as follows: 

1. Construct the N x R stoichiometric matrix, S, 
whose i — j element is given by — . 

2. Construct the macroscopic rate function vector / 

with elements fj = kj rim=i([^m]) Smj an( i tne di- 
agonal matrix F with elements Fa = f t . 

3. Construct the Jacobian matrix J whose i — j ele- 
ment is given by Vj(S ./),. Construct the diffusion 
matrix D = S ■ F • S T . 

4. The stochastic differential equations (linear 
Langevin equations) approximating the chemical 
master equation for the set of reactions (12) in the 
limit of large molecule numbers are given by 4 

j t m = j • m + ^- i/2 s ■ vTf W, (14) 

where 77,; (t), the iih entry of the vector rf(t), 
denotes the fluctuations about the macroscopic 
steady-state concentration of species i, i.e., r)i(t) = 
(rii(t)/Q) — [Xi\. The R dimensional vector T(t) 
is white Gaussian noise defined by (Ti(t)) = and 
(T i (t)r j (t')) = 6 i ,jS(t-t'). 

5. The covariance matrix g_ of the fluctuations in 
Eq. (14) is obtained by solving the Lyapunov 
equation 4 ' 5 

J • 3. + 3L ■ 1 T + D/fi = 0, (15) 

where ciy = (riiTij). The variance of the fluctua- 
tions is hence given by the diagonal elements of a . 
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C. Parameter values used in stochastic simulations 

Parameter values for the stochastic simulations shown 
in Fig. 1(a) are as follows. The enzymes, Malate dehydro- 
genase and Chymotrypsin, both have /3 = 1 ; we have sim- 
ulated only the first of these enzymes by using £7 = 17 fL, 
[E T ] = WnM, k 2 = 5s" 1 , k x = 5s" 1 , and k = 5 x 
10 7 M~ 1 s~ 1 . For the case P = 2.8 (the enzyme Lactate 
dehydrogenase), we used Q = 0.017 fL, [E T ] = 10 uM, 
k 2 = 210 s" 1 , k x = 75 s" 1 , and k = 3 x 10 6 M~ 1 s~ 1 . 
For the case j3 = 10 (a case compatible with a highly 
efficient enzyme), we used il = 170 fL, [Et] = InM, 
k 2 = Is- 1 , ki = OAs- 1 , and k = 10 8 M~ 1 s~ 1 . In all 
cases the total number of enzyme molecules tit was 100. 
Parameter values for Fig. 1(b) are f3 = 10, k 2 = Is" 1 , 



ki = 0.1s" 1 , [E T ] = InM, and fc = 10 8 M _1 s" 1 . 

The rate constants for the cases (3 = 1 and /3 = 2.8 
were obtained from the experimental studies. 6 " 8 The rate 
constants for j3 = 10 were not for a specific enzyme and 
hence were chosen from the known physiological ranges: 
for k 2 the range is 9 1 — 10 4 s" 1 , for fco the range is 10 
10 6 - 10 8 s" 1 M" 1 , and for K M the range is 9 10" 1 - 
10" 7 M. Similarly, the total enzyme concentrations were 
chosen from the physiological ranges: nano to millimolar 
concentrations. 11 The compartment volumes for the data 
in Fig. 1(a) were chosen such that the total number of 
enzyme molecules ut was 100 in all cases; for Fig. 1(b) 
the volumes were chosen such that nx could be varied 
over the range 1 to 100. 
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